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Abstract. A review of non-diffusive transport in fluids and 
plasmas is presented. In the fluid context, non-diffusive 
chaotic transport by Rossby waves in zonal flows is stud- 
ied following a Lagrangian approach. In the plasma physics 
context the problem of interest is test particle transport in 
pressure-gradient-driven plasma turbulence. In both systems 
the probability density function (PDF) of particle displace- 
ments is strongly non-Gaussian and the statistical moments 
exhibit super-diffusive anomalous scaling. Fractional dif- 
fusion models are proposed and tested in the quantitative 
description of the non-diffusive Lagrangian statistics of the 
fluid and plasma problems. Also, fractional diffusion opera- 
tors are used to construct non-local transport models exhibit- 
ing up-hill transport, multivalued flux-gradient relations, fast 
pulse propagation phenomena, and "tunneling" of perturba- 
tions across transport barriers. 



1 Introduction 

The widely used advection-diffusion equation rests on the va- 
lidity of the Fourier-Fick's prescription which in the case of 
transport of a single scalar, T, in a one-dimensional domain 
states that, q = —xdxT-\-VT, where q is the flux, x is the 
diffusivity, and V the advection velocity. From the statistical 
mechanics point of view, this model assumes an underlying 
Markovian, Gaussian, uncorrelated stochastic process. How- 
ever, despite the relative success of the diffusion model, there 
are cases in which this model fails to describe transport, and 
an alternative description must be used. The goal of this pa- 
per is to review some recent results on non-diffusive trans- 
port of particular interest to fluids and plasmas. We focus 
on non-diffusive Lagrangian particle transport and non-local 
transport of passive scalar fields. 



In the paradigmatic case of the Brownian random walk, the 
Gaussian statistics of the individual particle displacements, 
and the lack of correlations and memory effects (Markovian 
assumption), lead to a Gaussian PDF of the net particle dis- 
placement, and to the linear in time scaling for the mean, 
M ^t, and the variance, ~ t. Based on these scaling, 
the transport coefficients are defined as F = lim^^oo M(t)/t 
and X = l™^^oo {t) 1 1. The signatures of non-diffusive be- 
havior in Lagrangian particle transport include non-Gaussian 
PDFs of particle displacements and anomalous scaling of the 
moments of the form M and a'^ ^t^, with ^ 7^ 1 and/or 
7 7^ 1. When 7 > 1 (7 < 1) transport is super-diffusive (sub- 



diffusive), see for example |B ouchaud] ( | 1 99Q| ) . In either case, 
the diffusion model cannot be applied because the effective 
diffusivity x is either 00 or zero. 

The study of non-diffusive Lagrangian particle transport 
presented here focuses on two systems of interest to geo- 
physical fluid dynamics and plasma physics. In the geophys- 
ical context we consider transport in quasigeostrophic zonal 
flows. Quasigeostrophic flows are 2-dimensional, rapidly ro- 
tating flows in which there is a gradient in the Coriolis force. 
These flows are relevant in the study of mesoscale dynam- 



ics in the atmosphere and the oceans, see for example Ped- 



loskyl ([T987j). Zonal shear flows occur naturally in nature; 
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two well-known examples are the Gulf Stream and the po- 
lar night jet above Antarctica. Barotropic perturbations of 
these flows give rise to low frequency instabilities known 
as Rossby waves that have a key influence on the dynam- 
ics and transport. Following del-Castillo -Negrete and Mor^ 
rison| ( [1993| ); |del-Castillo-Negrete| ([1998j) we study chaotic 
transport by Rossby waves in zonal shear flows. In the 
plasma physics context we consider non-diffusive transport 
in pressure-gradient-driven plasma turbulence. This system 
is of relevance to the understanding of magnetically confined 
fusion plasmas. In this case, the Lagrangian particle dynam- 
ics corresponds to the motion of test particles in the presence 
of an external fixed magnetic field and a fluctuating turbulent 
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electrostatic potential. In the fluid and the plasma physics 
problems, we present numerical evidence of non-diffusive 
transport. In particular, in both cases, the PDFs of particle 
displacements are strongly non-Gaussian and the variance 
exhibits anomalous scaling of the super-diffusive type. 

As mentioned before, when the statistical moments ex- 
hibit anomalous scaling, the advection-diffusion model can 
not be applied and alternative models must be used. In this 
paper we review the use of fractional derivatives to construct 
such alternative models. Fractional derivatives are integro- 
differential operators that generalize the concept of deriva- 
tives of order n, to fractional orders [Samko et. al| ( |1993| ); 
|Podlubny| (|19991). Although the origins of fractional cal- 
culus go back to the origins of regular calculus, the use of 
fractional derivatives in the mathematical modeling of trans- 
port is relative recent. We present a brief review of this 
formalism in connection with the continuous time random 
walk (CTRW) model. The CTRW generalizes the Brown- 
ian random walk by incorporating non-Gaussian jump PDFs 
and non-Markovian waiting time PDFs Montroll and Weiss 
; [M'ontroU and Shlesinger|(i984| ); |Metzler and Klafter^ 
2000]). Following this, we construct effective macroscopic 
fractional diffusion models of the PDFs of particle displace- 
ments „del-Castillo-Negrete et.al| p004[ |2005| ). A compari- 
son is presented between the analytical solutions of the frac- 
tional models and the numerical results obtained from the 
Lagrangian statistics for the fluid and plasma problems men- 
tioned above. 

The use of fractional derivatives in transport modeling is 
close related to the problem of nonlocal transport. By nonlo- 
cal we mean that the flux of the transported scalar at a point 
depends on the gradient of the scalar throughout the entire 
domain. The generic mathematical structure of the nonlocal 
flux is q = —xj^{x — y)dyTdy, where the function JC mea- 
sures the degree of nonlocality. The "width" of this function 
depends on the strength of the non-locality, and in the limit 
when /C is a Dirac delta function, the flux reduces to the local 
Fourier-Fick's prescription. Motivated by the successful use 
of fractional derivatives to model non-diffusive Lagrangian 
transport, we discuss the use of these operators to construct 
non-local model of passive scalar transport. Following [del- 
Castillo-Negrete| ( |2006| ); |del-Castillo-Negrete, et. al| ( [2008j 



we present numerical results illustrating important non-local 
transport phenomenology including: up-hill transport, mul- 
tivalued flux-gradient relations, fast pulse propagation phe- 
nomena, and "tunneling" of perturbations across transport 
barriers. 

The rest of this paper is organized as follows. Section II 
discusses non-diffusive chaotic transport by Rossby waves in 
zonal flows. Non-diffusive turbulence transport in plasmas 
is studied in Sec. III. Section IV presents a brief review of 
fractional diffusion. The applications of fractional diffusion 
to model the PDFs of particle displacements in the Rossby 
waves and the plasma problems are discussed in Sec. V. Non- 
local transport is studied in Sec. VI, and Sec. VII presents the 



conclusions. 



2 Non-diffusive chaotic transport by Rossby waves in 
zonal flows 

In this section we study non-diffusive chaotic transport by 
Rossby waves in zonal shear flows. Since the flow is 2- 
dimensional and incompressible, the flow velocity can be 
written as v = (— Q^^,^^;^) where is the stream- 

function. In this case, the Lagrangian trajectories of individ- 
ual tracers, dr/dt = v, are obtained from the solution of the 
Hamiltonian system. 



dx _ dy _ 

dt dy dt dx 



(1) 



where ^ plays the role of the Hamiltonian and the r={x^y) 
spatial coordinates play the role of canonically conjugate 
phase space coordinates. Hamiltonian systems of the form 
in Eq. ([T]) are always integrable when ^ does not depend on 
time. However, when ^ depends explicitly on time, the sys- 
tem can be non-integrable and individual trajectories can be 
chaotic, see for example [Tabor ( 1989'). The main goal of the 
study of chaotic transport is to understand the global trans- 



port properties of tracers in this case, see for example pttino 
( |1989| ). Problems of particular interest to geophysical flows 
include the study of the formation and destruction of trans- 
port barriers |del-Castillo-Negrete and Morrison ([1993 ), and 
the study of the Lagrangian statistics del-Castillo-Negrete 
(1998 ). Here we focus on the second problem. 

To construct a model for the streamfunction ^f{x^y^t) we 
have to consider the dynamics of the system. In the case of 
quasigeostrophic flows, ^{x^y^t) is obtained from the poten- 
tial vorticity conservation law 



dt 



-(v-V)^ = 0, 



(2) 



where according to the /3 -plane approximation, = + 
py. We have adopted a right-handed Cartesian coordinate 
system with z pointing in the direction of the rotation of the 
system and y in the direction of the Coriolis force gradient. 
That is, y points in the "northward" direction and x is a peri- 
odic coordinate in the "eastward" direction. 

To simplify the solution of the non-linear Eq. ([2]) we as- 
sume a streamfunction of the form 



^ = ^Q{x,y)^^i{x,y,t), 



(3) 



where is the superposition of a zonal shear flow with 
dependence u^iy) = tanh?/, and a regular neutral mode in its 
co-moving reference frame, 

^0 = — In(coshy) +ei0i(?/)cos(/cix) +ci?/. (4) 
The function is a time dependent perturbation of the form 

^1 =e2(/>i(y)cos(/cix-a;t), (5) 
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where ei and 62 are free parameters determining the ampli- 
tude of the linear Rossby waves, and uj is the frequency of 
the perturbation. The eigenf unction ^ 



[1 +tanh?/] 



(l-ci)/2 



l — tdLuhy] 



(l+ci)/2 



(6) 



is obtained from the solution of the linear eigenvalue problem 
of the quasigeostrophic equation and (/ci,ci) are obtained 
from the corresponding dispersion relation for neutral (zero 
growth rate) modes del-Casti llo-Negrete| ( [ 1 998| ) . 

When 62 = the streamfunction is time independent and 
the solution of Eq. ([T]) can be reduce to a quadrature. In this 
case the Lagrangian dynamics is integrable and the orbits of 
the tracers can be classified in two types: (i) trapped orbits 
that encircle the vortices and (ii) untrapped orbits that move 
freely in the East- West direction following the zonal shear 
flow. These two types of orbits are separated by the separa- 
trix that joints the hyperbolic stagnation points of the flow. 
When 62 7^ 0, the system ceases to be integrable. In particu- 
lar, as shown in Fig. 1, the perturbation breaks the separatrix 
and creates a stochastic layer where tracers alternate chaot- 
ically between following the zonal flow and being trapped 
inside the vortices. 

To characterize transport in the chaotic regime we follow a 
statistical approach. The most basic quantity is the probabil- 
ity density function (PDF) of particle displacements. Trans- 
port in the "north- south" direction is trivial since particle or- 
bits in the ^/-direction are bounded by the zonal flows. There- 
fore, we focus on transport in the "east-west" direction, i.e. 
along the zonal flow. Given an ensemble of initial conditions, 
{{xi^yi)} with i = l,2,...A/p we compute the PDF of par- 
ticle displacements, P{6x^t) where Sxi{t) = Xi{t) — Xi{0). 
By definition, at t = the PDF is a Dirac delta function, 
P{Sx^t = 0) = S{Sx). As t increases, the PDF widens and 
might drift to one side or the other. Note that, although Sx is 
a periodic function in the annular domain shown in Fig.l, to 
compute the statistics we treat 6x as variable defined on the 
(—00,00) domain. 

To study the self- similar evolution of the PDF we intro- 
duce the scaling variable 

r] = {Sx-{Sx))t-^^\ (7) 

Figure 2 shows the rescaled PDF, t^^'^P, as function of 77. 
The observation that the rescaled PDFs collapse for succes- 
sive times leads support to the assumption that, at large times, 
P converges to a self-similar distribution of the form 

P*ix,t)=t-'</^fi7j), (8) 

where / is a scaling function, and 7 is the scaling exponent. 
The scaling in Eq. ^ implies the following scaling of the 
moments 

(X^)^t^7/2^^ (9) 

where X = Sx — {5x). Equation ^ also implies 

P*(X,t) = A^/^P* (A^/2X,At) , (10) 




Fig. 1. Chaotic transport by Rossby waves in the quasigeostrophic 
zonal flow in Eqs. ([TJ and ([3]). In the presence of two or more 
Rossby waves, the trajectories of passive tracers are typically 
chaotic. In particular, as shown in the figure, tracers alternate in 
a seemly random way between being trapped in the vortices and 
moving freely along the "east- west", x angular direction, following 
the shear flow flanking the vortices. 



where A is a free parameter. According to this relation, 
up to an scale factor, the limit distribution, P*, is invari- 
ant under the space-time renormalization operation (X,t) 
( A^/^X, Xt) . That is, the PDF at a later time can be obtained 
from a rescaling of the PDF at an early time. 

In the diffusive case, P* is a Gaussian, 7 = 1, and Eq. ([8| 
corresponds to the similarity solution of the advection- 
diffusion equation. However, in the numerical results shown 
in Fig. 2, transport is non-diffusive because 7 7^ 1 and the 
scaling function is not a Gaussian. In particular, the tails of 
the PDFs exhibit a decay significantly slower than Gaussian 
and a strong asymmetry. Because, 7 > 1, it is concluded that 
"East- West", azimuthal chaotic transport by Rossby waves 
is zonal flows is super-diffusive. For further details on the 
statistics and a dynamical explanation of the dependence of 
the asymmetry of the PDF on the perturbation frequency u 
see Ref. del- Castillo-Negrete| ([1998). This reference also 
discusses the comparison of the model presented here with 
experimental results on transport in rapidly rotating fluids 
.Solomon et. aU1993j . 

It is interesting to mention that there is a very close anal- 
ogy between the dynamics of Rossby waves in rapidly ro- 
tating neutral fluids in the quasigeostrophic approximation 
and drift- waves in magnetized plasmas, see for example 
[Petviashvili and Pokhotelov (19921); |Horton and Hasegawa] 
^1990j ; ,Horton and Ichikawa,(1996| ). In this analogy, the role 
of the rapid rotation is played by the strong magnetic field, 
the fluid streamfunction corresponds to the electrostatic po- 
tential, the fluid vorticity to the plasma density, and the gra- 
dient in the Coriolis force corresponds to the plasma density 



4 



D. del-Castillo-Negrete: Non-diffusive, non-local transport in fluids and plasmas 





Fig. 3. Fluctuating electrostatic potential 5> at a fixed time ob- 
tained from the numerical integration of the plasma turbulence 
model in Eqs. (T2])-|T4|. The dark (light) coherent patches de- 
note rotating (counter rotating) E x B eddies. The trapping ef- 
fects of these eddies along with intermittent large radial displace- 
ments caused by avalanche-like plasma relaxation events, give rise 
to non-diffusive transport and to the non-Gaussian PDF in Fig. 4 
|del-Castillo-Negrete et.al| ( [20041 [2005] ). 



Fig. 2. Rescaled probability distribution function (PDF), t^^^P, of 
passive tracers displacements, 6x{t) = x(t) — x(t = 0), as function 
of the similarity variable, 77 = {Sx — (Sx)) jV^^ with 7 = 1.9. The 
dynamics corresponds to the quasigeostrophic model in Eqs. ([TJ and 
The plot shows the PDF at t = 800, 900 and 1000. Consis- 
tent with the self-similar scaling in Eq. ([8]), the PDFs at successive 
times collapse. The anomalously large displacements induced by 
the zonal flow (see Fig. 1) result in the strong departure of the 77 < 
tail from the Gaussian fit (dashed line). The value 7 > 1 indicates 
super-diffusive transport. 



background gradient. Based on this analogy, as discussed in 
|del-Castillo-Negrete~[ ( |2Q00| ), the results presented here are 
directly applicable to the study of non-diffusive chaotic trans- 
port by drift waves in magnetized plasmas. 



3 Non-diffusive turbulent transport in plasmas 

In the example discussed in the previous section, transport 
resulted from chaotic advection. That is, from the chaotic 
dynamics of the deterministic equations describing the parti- 
cle orbits. In particular, the streamf unction ^ is a determin- 
istic differentiable function. In the case of turbulent transport 
the situation is different since the flow velocity advecting the 
tracers is a nondeterministic, random function. Nevertheless, 
turbulent systems can also exhibit non-diffusive transport of 
passive tracers. In this section we present an example in the 
context of plasma physics. 

As in the previous section, we follow a Lagrangian ap- 
proach and consider the statistics of a large ensemble of 



tracer particles. In the plasma, the particle motion responds 
to the combined effect of a turbulent electric field, E = 
— V<l>, and a fixed external magnetic field, Bq. The equation 
of motion of the tracers are obtained from Newton's law with 
the Lorentz force. However, in the guiding center approxi- 
mation, see for example [Nicholson| ( |1983| ), the equations can 
be simplified as the first order system 



2V^xBo, 



(11) 



where r = (x,?/) denotes the position of the particle in the 2- 
dimensional plane perpendicular to the magnetic field. This 
system has also a Hamiltonian structure with the potential, 
i>, playing the role of Hamiltonian. 

The fluctuating plasma electrostatic potential, <l>, is ob- 
tained from the solution of the turbulence model. Here, fol- 



lowing jdel-Castillo-Negret^ we consider 

pressure-gradient-driven turbulence in cylindrical geometry. 
The underlying instability of this type of turbulence is the 
resistive interchange mode, driven by the pressure gradient. 
This instability is the analogue of the Rayleigh-Taylor insta- 
bility responsible for the gravity-driven overturning of high 
density fluid laying above a low density fluid. In magneti- 
cally confined plasmas, the role of gravity is played by the 
magnetic field lines curvature. The turbulence model |Car-| 



reras, et.al ( 1987 ) is based on an electrostatic approximation 
of the reduced resistive magneto hydrodynamic equations. 
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Fig. 4. Rescaled probability distribution functions (PDFs), t^P, of 
passive tracers radial displacements, x(t), as function of the sim- 
ilarity variable, xjt^ with z/ = 2/3. The dynamics correspond to 
the pressure-gradient-driven plasma turbulence model in Eqs. |T2| - 
([T4|. The plot shows the PDF at t = 0.2, 0.4, 0.6 and 0.88. Like 
in the fluid dynamic case in Fig. 2, the collapse of the PDFs at suc- 
cessive times indicates a self-similar scaling of the form in Eq. ([8]). 
In this case, the non-diffusive transport manifest in the slowly de- 
caying non-Gaussian tails of the PDF. The value > 1/2 indicates 
super-diffusive transport fdel-Castillo-Negrete et.al| ( [20041 [2005 ) ). 



dt 



dr r dO 



(13) 



where <i> is the electrostatic potential, p the pressure, and 
d/dt = dr-\-^-S/. The instability drive is the flux- surface 
averaged pressure gradient, d{p)/dr, determined according 
to 



d{p) ^ 1 d ^ 
dr r dr 



Vrp)=So + D- 



1 d ( d{p)\ 



r dr \ dr J 



(14) 



The tildes indicate fluctuating quantities (in space and time), 
and the angular brackets, (), denote flux surface averaging 
over a cylinder at a fixed radius. The equilibrium density is 
no, the ion mass is m^, the averaged radius of curvature of 
the magnetic field lines is Tc, and the resistivity is rj. The sub 
indices "±" and "||" denote the direction perpendicular and 
parallel to the magnetic field respectively. The function 5*0 
represents a source of particles and heat which we model us- 
ing a parabolic profile, = [l — (r/a)'^] . Figure 3 shows 
a snapshot in time of the fluctuating electrostatic potential 4> 



obtained form the solution of Eqs. (^2])-(p4|). 

Having computed <l>, the next step is to integrate Eq. ( jj_ 
to obtain the orbits of the tracers. The initial condition con- 
sists of 25 X 10^ particles with random initial positions in 



and z, and radial position r = 0.5a. By definition, at 
t = 0, the PDF, P, of radial particle displacements, x = 
[r{t) — r(0)]/a, is a Dirac delta function. As time advances 
the P(x,t), spreads and develop slowly decaying, "fat" tails. 
Figure 4 shows the long-time time behavior of the PDF as 
function of the similarity variable x/t^. The strong non- 
Gaussianity of P is evident. Like in the previous fluid exam- 
ple case, transport is super-diffusive because v > 1/2. Ev- 
idence of non-diffusive transport has also been observed in 
other plasma systems including gyrokinetic simulations of 



ion temperature gradient (ITG) turbulence Sanchez, et. al 
C2008j ). 

4 Fractional diffusion models of non-diffusive transport 

One of the main goals of transport modeling is to construct 
effective macroscopic transport equations that reproduce ex- 
perimentally or numerically observed phenomena. For ex- 
ample, in the fluid and the plasma transport problems dis- 
cussed in the previous two sections, the goal is to con- 
struct a transport equation that describes the observed spatio- 
temporal evolution of the PDF, P, of particle displacements. 

When transport is diffusive, a simple solution to this prob- 
lem is provided by the advection-diffusion equation 



(15) 



where the advection velocity and diffusivity are obtained 
from the asymptotic behavior the statistical moments 



V = lim 



t 



X-- 



lim 



{x^ 
2t 



(16) 



of the particle's displacements, x. However, this approach 
fails in the case of non-diffusive transport. In particular, 
according to the scaling in Eq. ^ when there is super- 
diffusion, X ^ cxo. Moreover, as it is well-known, the 



Green's function of Eq. ( 15 ) in an unbounded domain, is a 
translated Gaussian and this significantly limits the range of 
PDFs that this model can describe. In particular, PDFs with 
7^1 scaling and/or with slowly decaying tails, like those 
obtained in the examples discussed before (Figs. 2 and 4), 
cannot be modeled using a simple advection-diffusion equa- 
tion. 

From the statistical mechanics point of view, the 
advection-diffusion model assumes an underlying Marko- 
vian, Gaussian stochastic process with a drift, i.e. a bi- 



ased Brownian random walk, see for example Paul and 



Baschnagel ( 1999). However, the description of transport in 
the presence of coherent structures requires the use of ran- 
dom walk models that incorporate more general stochastic 
processes. In particular, in the fluid problem discussed in 
Sec. 2, the trapping effect of the vortices gives rise to non- 
Markovian effects, and the zonal shear flows give rise to non- 
Gaussian particle displacements. In the plasma physics prob- 
lem discussed in Sec. 3, the non-Markovian effects are due 
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to the trapping in electrostatic eddied, and the non-Gaussian 
particle displacements result from avalanche-like radial re- 
laxation events. 

The Continuous Time Random Walk (CTRW) model 
MontroU and We iss (196^; 



MontroU and Shlesinger| ( p^84l ); = iie 



Metzler and Klafter| ( |20QQ| ) provides an elegant powerful 
framework to incorporate these type of effects. The CTRW 
generalizes the Brownian walk in two ways. First, contrary 
to the Brownian random walk where particles are assumed 
to jump at discrete fixed time intervals, the CTRW model 
allows the possibility of incorporating a waiting time prob- 
ability distribution, In addition, the CTRW model al- 
lows the possibility of using non-Gaussian jump distribution 
functions, 77 (x), with divergent moments to account for long 
displacements known as Levy flights. Given ip and 77, the 
probability of finding a tracer at position x and time t is de- 
termined by the MontroU- Weiss master equation 



dtP = 



/ dt'(j){t-t') / dx'[ri{x-x')P{x',t')- 

Jo J -co 



(17) 



-r]{x-x')P{x,t')], 



The spatial integral on the right-hand- side represents the 
gain-loss balance for P at x. In particular, the first term in- 
side the square bracket gives the increase of P due to par- 
ticles moving to x while the second term describes the de- 
crease of P due to particles moving away from x. The time 
integral accounts for memory effects weighted by the func- 
tion (j){t). In Fourier-Laplace variables. 



Akx 



r]{x)dx, 



C[(l)]=4>{s)= / e''(l){t)dt, 
Jo 



Eq. ( 17 ) takes the form 



P{k,s). 



1 



1 —ip{s)fi{k) 



(18) 



(19) 



(20) 



where the relation between the waiting time PDF and the 
memory function is ^ = sij) / (^1 — ijj ^ . 

The MontroU-Weiss master Eq. (Tl) can be used directly 
to model non-diffusive transport, see for example |van Mil-| 
ligen et. al| ( |2Q04| ); [Spizzo et. al| ( [2009 ). However, this de- 
scription carries in a sense too much information concerning 
the details of the underlying stochastic process that might be 
irrelevant in the long-time, large-scale description of trans- 
port. This motivates the derivation of a macroscopic trans- 



port equation from Eq. (20) valid in the time asymptotic 
(s 0) long-wavelength (k 0) "continuum" limitjSaichev 
and Zaslavsky| ( pW7l ); [Metzler and Klafter| ( [2000| ); |Scalas et. 



al|( |2004| ). A key aspect of this limit is that only the asymp- 
totic behavior, i.e., the tails of the r] and ijj PDFs matter. This 
is a significant advantage over the use of the kinetic master 



equation that requires the detailed knowledge of these func- 
tions. 

As expected, in the Markovian- Gaussian case 

1 



-fit 



r]{x) - 



271(7 



(21) 



where (t) = is the characteristic waiting time and = 
(x^) is the characteristic mean square jump, the fluid limit of 



the master equation Eq. ( 20 ) leads to the standard diffusion 
equation in ( p?] ). However, the situation is quite different in 
the case of algebraic decaying PDFs of the form 



Tj ^ \X\ 



-(a+1) 



(22) 



where for simplicity we have assumed that r] is symmetric. 
In this case, for < ^5 < 1, {t) diverges, and there is no char- 
acteristic waiting time. Similarity, for a < 2, (x^) diverges, 
indicating a lack of characteristic transport scale. The use 
of this type of algebraic decaying PDFs is motivated by the 
significant probability of very large trapping events and very 
large spatial displacements, as it is the case in the examples 
discussed in Sees. 2 and 3. From the asymptotic behavior in 
Eq. ([22|) it follows that for small s and k. 



7}(fe)^i-ifer+. 



(23) 



Substituting Eq. ( 23 ) into Eq. d20| we get to leading order 



s^P{k,s) 



-xWP^Ks). 



(24) 



To obtain the macroscopic transport equation we need to in- 



vert the Fourier-Laplace transforms in Eq. (24). This can be 
formally done by writing 

IDP,P = xD%^P, (25) 
where the operators in Eq. ([25]) are defined according to 



C ID^P =s^P{x,s)-s^-'8{x), 



WP{k,t), 



(26) 



(27) 



for < /3 < 1. Equations ( [26] ) and ( [27] ) are the natural gener- 
alizations of the Laplace transform of a time derivative and 
the Fourier transform of a spatial derivative. This motivates 
the formal identification of the operator qD^ as a "fractional 
time derivative" for < /3 < 1, and the operator D^^^ as a 
"fractional space derivative" for 1 < a < 2. As expected, for 
aov P integers, the regular derivatives are recovered. 

The previous discussion assumed a symmetric jump 
stochastic process, iri{x) = r]{—x). It can be shown that in 
the general case the transport equation is 

gZ)fP = X[;-ocD^+r.Z)^]P, (28) 

where the operators on the right hand side are the left and 
right Riemann-Liouville fractional derivatives of order a 
Samko et. al[ ( [T993l ); [Podiub5^ ( [T^ 

1 5- r f{y) 



f - 



T{m — a) dx^ 



{x-y)^ 



+ l-m 



(29) 
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i-iy 



r(m- 



f 



fiy) 



iy-x) 



-dy, 



(30) 



where m is a positive integer such that m — l<a<m. In 
this general formulation, the asymmetry of the underlying 
stochastic process manifests on the parameters / and r, 



2cos(a7r/2) ' 



(1 + ^) 
2cos(a7r/2) ' 



(31) 



that control the relative weight of the left and right fractional 
derivatives, where — 1 < ^ < 1. In the symmetric case, ^ = 0, 
+ xD'^] which corresponds to the 



-^\x\ ~ 2cos(7ra/2) 



operator defined in Fourier space in Eq. ( [27] ). In the time 
domain, the fractional derivative operator in time, o^f , in- 



troduced in Eq. ( 26 ) become an integro-differential operator 
of the form 



1 



T{l-P)Jo (t-t'Y 



dt' , 



(32) 



where < /3 < 1. For a derivation of fractional diffusion 
models that incorporate more general stochastic processes, 
including the physically important case of truncated Levy 
statistics, see |Cartea and del-Castillo-Negrete| (j2007|. For 
a derivation of fractional diffusion models using quasi-linear 



type renormalization techniques see Sanchez, et. al (2006). 



5 Applications of Fractional diffusion models 

The goal of this section is to use the fractional diffusion equa- 
tion to model the non-diffusive transport of tracers discussed 
in Sees. 2 and 3. In particular, we show that the numerically 
obtained PDFs of the particle displacements in Figs. 2 and 
4 can be obtained as solutions of effective macroscopic frac- 
tional diffusion equations. 



The solution of the initial value problem of Eq. (28 ) with 

P(x,t = 0)=Po(^) is 



P{x,t)-- 



Po{x)G{x-x,t)dx , 



(33) 



where the Green's function (propagator) G is the solution of 
the initial value problem G{x^t = 0) = S{x) with S{x) the 
Dirac delta function. Using Eqs. ([26]) and ( [27] ), the Fourier- 
Laplace transform of Eq. ([28]) leads to the solution 



G = 



where 



(34) 



(35) 



for a / 1. Introducing the Mittag-Leffler function, see for 
example [Podlubnyj ([1999]) , 



n=0 



r(/3n + l)' 



C[E^{ct^)] 



(36) 



the inversion of the Fourier-Laplace transform in Eq. (34) 
gives 

Gix,t)=t-^/''Kirj), 



where 

7? = x(x'/^t)-^/^ 



(37) 
(38) 

(39) 



is the similarity variable. Further details of the solution of 
the initial value problem and useful asymptotic and con- 
vergent expansions of the Green's function can be found 
in Refs. [Metzler and Klafter] ( [20001 ); [Saichev and Zaslavsky| 
( [T997] ); [Marnardi et.al| ( [20QT] ). 

Of particular interest is the asymptotic behavior in x, for a 
fixed t = to, 

G{x,to)^x-^^+^\ x» (xy^to)^^". (40) 
and the small t and large t scaling at fixed x = xo, 



G{xo,t)^ < 



for t<C (x/^^o) 
t-^for t>(^x7^^o) 



(41) 



From these relations it follows that the order of the fractional 
derivative in space, a, determines the algebraic asymptotic 
scaling of the propagator in space for a fixed time, and the 
order of the fractional derivative in time, (3, determines the 
asymptotic algebraic scaling of the propagator in time for a 
fixed X. These two properties provide a useful guide to con- 
struct fractional models given the spatio-temporal asymptotic 
scaling properties of the PDF. Using Eq. ( [37] ), the moments 
in the fractional model are given by 



that implies the anomalous diffusion scaling 



(42) 



(43) 



According to Fig. 2, the scaling exponent of the PDF of 
particle displacements in chaotic transport by Rossby waves 
is 7 ~ 1.9. As expected, this value is also consistent with 
the scaling of the second moment computed directly form 
the Lagrangian statistic of displacements. Based on this, in 
the construction of the fractional model we assume 7 = 2, 
which according to Eq. ( [43] ) implies a = (3. This special case 
corresponds to the neutral fractional diffusion equation, for 
which G in Eq. ([37]) is |Mainardi etaT| ( |2001J : 



G{x,t)-- 



r 



sin[7r(a-C)/2]<- 



TT l + 27y"cos[7r(a-C)/2]+7?' 



2a 



(44) 
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the PDFs of particle displacements according to Eqs. ([40) 



-800 -600 -400 -200 

Sx 



200 



Fig. 5. Comparison between the PDF of particle displacements, Sx, 
in the quasigeostrophic zonal with Rossby waves(solid line), and the 
PDF obtained from the solution of the fractional diffusion model in 
Eq. (|28j with a = /3 = 0.9, and 0=1. 




Fig. 6. Comparison between the PDF of particle displacements, x, 
in the resistive, pressure-gradient-driven plasma turbulence model 
in Eqs. |T2])-|T4| and Fig. 4 (triangles), and the PDF obtained from 
the solution of the fractional diffu sion model in Eq. ([28|) wit h a = 
3/4,1 3 = 1/2, 6* = 0, and x = 0.09 |del-Castillo-Negrete et!al|p004 
[2005). 



for 77 > where r] = Sx/t is the similarity variable and 
6 = tan(7rC/2)/tan(7ra/2). The solution for 7^ < is ob- 
tained using the relation K{—r];aX) = K{r];a^—Q. Fig- 
ure 5 shows a comparison between the fractional diffusion 



solution in Eqs. ( 33 ) and (44 ) with the PDF obtained in Sec. 2 
from the Lagrangian statistics of the quasigeostrophic trans- 
port problem. 

In the case of turbulent transport in pressure-gradient- 
driven plasma turbulence, the asymptotic scaling analysis of 



and (41 ) and the super-diffusive scaling of the moments in 
Eq. (42), indicate that a = 3/4 and /3 = 1/2. Figure 6 com- 



pares the solution of the fractional diffusion equation for 
these parameters with the PDF obtained from the direct nu- 
merical simulation shown in Fig. 4. Details on the explicit 
solution of the fractional diffusion equation can be found 
del-Castillo-Negrete et.aT[ (|2QQ4| [2005). As discussed in 



m 
Sec 



3, the Lagrangian study of transport in plasmas was 
based on the guiding-center equations of motion which are 
an approximation to the dynamics valid in the limit of zero 
Larmor radius. The role of finite Larmor radius effects on 
non-diffusive transport, an in particular on fractional diffu- 
sion was studied in Ref . jGustaf son et. aT| ( |2QQ8| ). 

6 Non-local transport 

In the previous sections we discussed non-diffusive transport 
in the context of test particle Lagrangian transport in fluids 
and plasmas. One of the main goals was to construct macro- 
scopic effective transport models to describe the PDF of par- 
ticle displacements in chaotic and turbulent flows. It was 
shown that fractional diffusion operators provide a frame- 
work to describe the spatio-temporal evolution of the PDFs. 
In particular, the long tails of the PDFs as well as the non- 
Gaussian scaling of the Lagrangian statistics are well capture 
by fractional diffusion models. Motivated by these results, in 
this section we discuss the use of fractional diffusion mod- 
els to describe non-diffusive transport of passive scalars, like 
temperature, density, pressure or the concentration of a pol- 
lutant in flow. 

The starting point is the conservation law 



dtT = -d,q, 



(45) 



where T denotes the scalar field transported and q denotes 
the flux. For simplicity we limit attention to the transport of 
a single scalar in a 1 -dimensional domain. The conservation 



law ( 45 ) has to be complemented with a prescription relating 
q and T. In the case of diffusive-transport this closure is 
provided by the Fourier-Fick's local prescription 



-xd^T^VT, 



(46) 



where x is diffusion coefficient and V is the advection veloc- 



ity. Substituting Eq. ( [46|) in to Eq. ( [45] ) leads to the advection- 
diffusion model in Eq. \\5\. 

Although the advection-diffusion model has been success- 
fully applied to a wide variety of transport problems, there 
are cases in which this model fails to describe the dynam- 
ics. The examples discussed before showed clear evidence 
of this in the case of the PDF of Lagrangian particle dis- 
placements. Here we explore the role of non-diffusive trans- 
port of scalars, like temperature, for which a Lagrangian test 
particle perspective might not be readily available. One of 
the main motivations for this study is the understanding of 
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fast propagation phenomena in magnetically confined plas- 
mas. The basic problem can be understood without entering 
into the details concerning the plasma system. The top panel 
in Fig. 8 shows the basic configuration of interest, where 
To is perturbed by a pulse at the edge of the domain. The 
problem then is to study the relaxation of the system back to 
the steady state. These type of perturbative transport experi- 
ments are commonly performed in magnetically confined fu- 
sion devices where a plasma is suddenly cooled at the edge. 
It has been observed in several experiments that such cold 
pulse perturbations travel from the edge to the center of the 
device at speeds significantly greater than the typical diffu- 
sive time scales. Because of this, attempts to model some of 
these experiments using the diffusion equation have failed. 
Here we discuss the use of non-local transport models as an 
alternative to diffusive models to describe these phenomena. 
By non-local we mean that, contrary to the Fourier-Fick's 



local prescription in Eq. (46), the flux g at a given point de- 
pends on the gradient of T throughout the entire domain. The 
generic mathematical structure of these nonlocal models is 



(47) 



<?(^) = -X / JC{x-y)dyT{y)dy, 



where the kernel /C determines the level of non-locality. In 
the case when JC = 5{x — x'), Eq. ( [47|) reduces to the familiar 
Fourier-Fick prescription in Eq. |46| ), where for simplicity 
we assume V = 0. 

Non-local transport is a problem of significant interest in 
plasma physics, see for example [Callen and Kissick] ( |1997| ) 
and references therein. Fux-gradient relations of the form in 



Eq. ( 47 ) have been used in the study of parallel electron heat 
transport in magnetized plasmas ( |Held et.al| ( |20Ql )), and in 
the study of transport due to long scale-length fluctuations 



( Yoshizawa et. al| ( |2QQ3| )). However, the physics behind the 
non-local models discussed here is different, and it is based 
on the theory of non-Gaussian stochastic processes. Moti- 
vated by the results discussed in the previous sections, we 



model the non-local flux-gradient relation in Eq. ( 47 ) using 
fractional derivative type operators of the form 



q = -x{x)[laV^-r,V^]T, 
where x can depend on x and 
1 



r(2-a)y„ (x-y) 



^^^^ "r(2-a) 



' r{b)-r {y) 

{y-x) 



a-l 



dy, 



where V 
that the operators 



dxT, and / and r are defined in Eqs. (31 



(48) 



(49) 



(50) 



Note 



and x'^b exactly the usual 



Riemann-Liouville fractional derivative operators introduced 
in Eqs. ( [29] ) and ([30|). As discussed in RefS|del-Castillo- 



ta 
)6 



Negrete] ( |2Q06l ); |del-CastiUo-Negrete, etTall ( [2Q08J this dif- 
ference has to do with the important issue that in a finite 



size domain x e (a^b) the Riemann-Liouville operators must 
be regularized to incorporate general boundary conditions of 
physical interest. 

The role of non-locality and asymmetry in transport is il- 
lustrated Fig. 7 that shows the time evolution of a localized 



pulse initial condition in the model in Eqs. ( [45] ) and ( [48] ) with 
a = 1.3 and 6> = 0.5. As the top panel shows, due to the 
asymmetry, 6^0, the peak of the distribution shifts to the 
right. It can be shown that the peak of the profile, Xm, during 



the drift satisfies [del-Castillo-Negretel ( [2006 ). 



where 



Vrr 



1 



2a 



(51) 



(52) 



As expected, in the 6> = symmetric case and in the a = 2 
diffusive limit, the drift vanishes. This drift results from 
the existence of "up-hill" transport which is a generic fea- 
ture of non-local transport models. In the Fourier-Fick's pre- 
scription the flux dependence on the local gradient is always 
"down-hill", i.e., in the direction opposite to the local gra- 
dient. However, as the vertical lines in the top and middle 
panels of Fig. 7 indicate, in this case there is region of "up- 
hill" transport in which the flux is in the same direction as 
the gradient. Moreover, as the bottom panel of Fig. 7 shows, 
in the non-local decay of the pulse, the flux-gradient rela- 
tion is not linear like in the Fourier-Fick's diffusive case, it is 
in fact multivalued. The multivalued relation between q and 
—dxT isa. generic feature of non-local transport models with 
or without asymmetry. 

In the study of the propagation of pulse perturbations, the 
first step is the computation of the steady equilibrium temper- 
ature profile, To{x) in the presence of a source of the form 



S = So exp 



{x-/j.sy 



2a? 



(53) 



with /ig = 0, and as = 0.075. For each simulation, the source 
amplitude was selected so that To(0) = 1. The simulations 
followed the spatio-temporal evolution of the perturbed tem- 
perature, ST{x^t) = T{x^t) — To(x), with initial condition 



ST{x,0) = -Aexp 



(54) 



where A = 0.3, /Hp = 0.75, and cr^ = 0.03. Details on the 
numerical method used to integrate the fractional transport 



model can be found in |del-Castillo-Negrete| ( [2006] ) . The bot- 
tom panel in Fig. 8 shows the time evolution of the normal- 
ized tracer perturbation, ST = ^T/|min[^T(x,0)] |, at differ- 
ent locations along the x-domain. We define the mean pulse 
propagation speed as the ratio of the normalized distance 
and the time delay, Vp = l/5t. The time delay is defined 
as the time required for the scalar field at x = to exhibit 
a drop of size STc. That is, ST{O^St) = STc. For the value 



10 



D. del-Castillo-Negrete: Non-diffusive, non-local transport in fluids and plasmas 



10 



T 5 



10 



(a) 











, 1 

1 






1 






1 




1 


1 

/J 




1 


// 




#^ 




i \ 
I \\ 









0.5 
X 



(b) 



Right 

Total 






^ 








Left 



0.5 
X 




non-locality can lead to a considerable increase of the pulse 
speed. In particular, the numerical results show that for the 
same value of Vp for a = 1.25 is about 10 bigger than 



the diffusive speed. This idea was used in Ref. del-Castillo 



Negrete, et. al| ( |2008| ) to model perturbative experiments 



on cold temperature pulse propagation in the Join European 
Torus (JET) magnetically confined controlled fusion device. 

To conclude we present recent results on the role of non- 
locality in the propagation of pulses through transport barri- 
ers. The local and non-local diffusivities are assumed to be 
of the form 



and 

Xnl '- 



XnlO 



tann I - 



L J 



)+tanh(|) 



(55) 



(56) 



-{x-xo) /w 



The tanh profile in Xni is introduced to guarantee the vanish- 
ing of the non-local flux in the core region where transport is 
assumed to be dominated by diffusive processes. The trans- 
port barrier is modeled by introducing a dip, e~^^~^^^ 
in the diffusivity profiles. In the calculations reported here 
Xdo = hxo = 0.5, C = 0.95, Xnio = hxc = 0.1, L = 0.025, 
and w = 0.005. In the non-local simulations, a = 1.25. Fig- 
ure 9 shows the spatio-temporal evolution of ST. The top 
panel shows the case of diffusive transport, Xnio = 0, in the 
absence of transport barriers. In this case, the pulse spreads 
throughout the plasma domain in a slow, diffusive time scale. 
As expected, as shown in the middle panel, in the presence 
of a transport barrier the diffusive propagation of the pulse 
is stopped. However, in the presence of non-local transport 
the pulse dynamics is fundamentally different. As the bot- 
tom panel in Fig. 9 shows, in this case the pulse can in fact 
go through the transport barrier. This "tunneling" effect is a 
unique novel property of non-local transport. 



Fig. 7. Non-local transport of a localized pulse initial condition ac- 
cording to the fractional diffusion model in Eqs. ( [45| ) and ( [48| l with 
a = 1.3, /3 = 1, and = 0.5. The solid line in panel (a) shows the 
profile at the final time, and the dashed lines the profiles at earlier 
times. The drift of the distribution results from the asymmetry ^ / 
of the fractional operator. Panel (b) shows the left, qi, the right, qr, 
and the total non-local flux, q, and panel (c) shows the flux-gradient 
relation. Contrary to the Fourier-Fick's linear relation, q — —xdxT, 
q and —dxT exhibit a nonlinear, multivalued relation. The top, left 
quadrant, —dxT<0 and q>0 corresponds to up-hill transport that 
occurs in the region bounded by the two vertical lines in panels (a) 
and (b). 



of the threshold we choose 5Tc = —0.0375. We considered 
three case: an a = 2 diffusive case, and two fractional cases 
with a = 1.7b and a = 1.25. The main conclusion is that 



7 Conclusions 

We have presented a review of recent results on non-diffusive 
transport in fluids and plasmas. The approach was based 
on the study of the Lagrangian statistics of large ensembles 
of particles. In general, the stochasticity in the Lagrangian 
trajectories can result from deterministic chaos or from tur- 
bulence. The examples discussed encompass both possibil- 
ities. In the studied of transport by Rossby waves in quasi- 
geostrophic zonal flows, the advection velocity was a smooth 
deterministic function but the Lagrangian trajectories exhib- 
ited Hamiltonian chaos. On the other hand, in the E x B 
transport plasma problem, the advection velocity was a non- 
deterministic random function obtained from the solution of 
a turbulence model. The main object of study was the proba- 
bility density function (PDF) of individual particle displace- 
ments, also know as the propagator. Both, the fluid chaotic 
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Fig. 8. Non-local fast pulse propagation. As shown in the top panel, 
perturbative transport simulations follow the evolution of a local- 
ized perturbation (dashed line) of an steady state passive tracer pro- 
file (solid line). The bottom panel shows the time traces of the nor- 
malized tracer perturbation, 5T — ^T/|min[^T(x,0)] |, at different 
locations along the x domain. In the local diffusive case (dashed 
line) the normalized propagation speed from the edge, x = 0.75, to 
the center, x = 0, of the domain is = 1. In the fractional case 
with a = 1.75 (solid line), Vp = 6.3, and in the fractional case with 
a = 1.25 (dotted line), Vp = 9.6. 



transport problem and the plasma turbulent transport prob- 
lem, exhibited strongly non-Gaussian spatio-temporal self- 
similar PDFs. In addition, the Lagrangian statistics in both 
cases exhibited super-diffusive scaling, <x^ >r^t'^ with 7 > 
1. The modeling of these PDFs using advection-diffusion 
equations is out of the question because the effective diffu- 
sivity diverges, and the propagators have non-Gaussian de- 
caying tails. The observed non-Gaussian statistics in the ex- 
amples discussed has its origin on the combination of anoma- 
lously large particle displacements, known as "Levy flights", 
and the trapping effects of coherent structures like fluid vor- 
tices and E X B plasma eddies. 



(a) 




X 



Fig. 9. Non-local "tunneling" of perturbations across a transport 
barrier. The figure shows the space-time evolution of the nor- 
malized passive tracer perturbation 5T = ^T/|min[^T(x,0)] | with 
dark blue (red) denoting 5T —1 {5T — 0). The top panel corre- 
sponds to diffusive transport in the absence of transport barriers. 
The middle and bottom panels correspond to diffusive and non-local 
transport respectively in the presence of a transport barrier The ver- 
tical dashed line indicates the location of the transport barrier. 



We have shown that the PDFs of particle displacements 
can be modeled using fractional diffusion equations in which 
regular derivatives are replaced by fractional derivatives. 
Fractional derivatives are integro-differential operators that 
provide a powerful, elegant framework to incorporate non- 
Gaussian and non-Markovian effects on transport models. 
These operators naturally appear in the continuum limit of 
generalized random walk models that extend the Brownian 
motion by allowing non-Gaussian jump distribution func- 
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tions and general waiting time distribution functions. 

Going beyond the study of non-Gaussian Lagrangian 
statistics, we discussed the application of fractional deriva- 
tives to model non-local transport. The cornerstone of the 
diffusive transport paradigm is the Fourier-Fick's prescrip- 
tion according to which the flux at a given point depends only 
of the gradient of the transported field at that point. On the 
other hand, in the case of non-local transport, the flux can 
depend on the gradient throughout the entire domain. Al- 
though in many cases transport problems follow the Fourier- 
Fick's prescription, there are important situations in which 
this is not the case. A clear example is the fast propagation 
phenomena observed in perturbative transport experiments 
in magnetically confined plasma fusion devices. Motivated 
by the successful use of fractional derivatives in the study of 
non-diffusive Lagrangian transport, we used these operators 
to construct non-local models of passive scalar transport. We 
presented numerical results illustrating important non-local 
transport phenomenology including: up-hill transport, mul- 
tivalued flux-gradient relations, fast pulse propagation phe- 
nomena, and "tunneling" of perturbations across transport 
barriers. 

Some of the results presented here pertain specific sys- 
tems, i.e., Rossby waves in zonal flows and pressure- 
gradient-driven plasma turbulence. However, it is important 
to realize that the observed non-diffusive phenomenology de- 
pends on very general non-Gaussian statistical properties and 
not on specific details. In particular, other systems with co- 
herent structures and/or strong spatio-temporal correlations 
are likely to exhibit similar non-diffusive and non-local trans- 
port dynamics. 
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